INTEGRABILITY OF DIFFERENTIAL DIFFERENCE EQUATIONS WITH 

DISCRETE KINKS 



CHRISTIAN SCIMITERNA AND DECIO LEVI 

Abstract. In this article we discuss a series of models introduced by Barashenkov, Oxtoby and 
Pelinovsky to describe some discrete approximations to the <j> i theory which preserve travelling 
kink solutions. We show, by applying the multiple scale test that they have some integrability 
properties as they pass the A\ and A2 conditions. However they are not integrable as they fail 
the A3 conditions. 



1. Introduction 

In a recent article Barashenkov, Oxtoby and Pelinovsky [lj considered those exceptional 
discretizations of the (j) 4 theory in which the stationary kink can be centered anywhere between 
the lattice sites. The request that the kink be translational invariance generates three families 
of exceptional discretizations which provide new differential-difference equations which are not 
known to be integrable. The continuous (f> 4 model in itself is not completely integrable, i.e. it has 
not an infinite number of generalized symmetries, conservation laws and exact solutions, however 
it is an approximation of many C or S integrable nonlinear models like the Liouville equation 
or the sine-Gordon equation. So the analysis of their integrability is an interesting problem 
worthwhile to be investigated. To do so we apply the multiple scale analysis to the nonlinear 
discrete equations as this provides an integrability test and gives an integrability grading of the 
equations. 

In Sections and [2] and [3] we present a review of the models presented in [I] and of the results 
necessary to apply the multiple scale analysis to them. Then in Section H] we present the results 
of the calculations and in the last Section we give some concluding remarks. The details of the 
calculations are left to an Appendix. 

2. Exceptional models 

(j) 4 theory has been one of the most important nonlinear models for the description of statisti- 
cal mechanics and field theory systems (2j[19] . The discrete analogs of the cfi 4 kinks are solutions 
of equations of the form 

.. u n+i ~ 2n n + n n _i , w n 

(1) Un = ^2 ~2~ + Qn{U n -i,U n , U n+1 ), 

where u n = u{hn,t), with h a constant lattice spacing. The function Q n is chosen so as to 
give in the continuous limit, at first order in h, the (f> 4 potential Q = — \v? . These kinks have 
been used to describe charge-density waves in polymers and metals |21j . narrow domain walls 
in ferroelectrics [5], discommensurations in dielectric crystals [8j, and topological excitations in 
hydrogen-bonded chains [3l[10]. Physically, one of the most significant properties of domain 
walls and topological defects is their mobility [lj. Mathematically, discrete equations ([rj are 
known to admit stationary kink solutions [HUH]; however whether traveling discrete kinks exist 
remains an open question [9jlll p2ll24H26j . The discretization breaks the Lorentz invariance and 
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the existence of traveling discrete kinks becomes a nontrivial matter. Indeed, if the equation 
is non authonomous, the discretization even breaks the translation invariance of ([T]) due to the 
presence of a Peierls-Nabarro barrier, an additional periodic potential induced by discreteness. 
Miraculously, there are several exceptional discretizations which, while breaking the translation 
invariance of the equation, allow for the existence of translationally invariant kinks, that is, 
kinks centered at an arbitrary point between the sites. One such discretization was discovered 
by Speight and Ward using a Bogomolny-type energy-minimality argument [24, 25]. In [1] we 
have a systematic study of these exceptional cases based on the observation that the translational 
invariance of the kink implies the existence of an underlying one-dimensional map u n +i = F(u n ). 
A simple algorithm based on this observation gives three classes of exceptional discretizations: 



where in the first two equations the coefficients (/xi, H2, ^3), as presented in [T], are not indepen- 
dent. 

3. Multiple scale expansion of lattice equations 

For completeness we present here the basic ideas of the multiple scale expansion of lattice 
equations as presented in [15] in the case of expansions in terms of analytic functions. 

3.1. Expansion of real dispersive partial difference equations. 

3.1.1. From shifts to derivatives. We consider a function u n : Z — > M. depending on a discrete 
index n£Z and suppose that: 

(a) The dependence of u n on n is realized through the slow variable n\ = en € M, e € R, 
< £ <C 1, i.e. u n = u(ni); 

(b) n\ varies in a region of the integer axis such that u (rii) is therein analytical; 

(c) The radius of convergence of the Taylor series in n\ is wide enough to include as inner points 
all the points involved in the discrete equation. 

Under these hypotheses we can write the action of the shift operator T n , defined by T n u n — 
u n+ \ = u{n\ + e) in the variable n\, as 





(2b) Q 2 




(3) 



T n u(ni) 



u 



(ni) + eu (1) (ni) + — u (2) (m) + ••• + -ru^Hm) + ••• 
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where u^(ni) = d = d l ni u(ni), being d ni the total derivative operator. From ([3]) it follows 
that we can write the following formal expansion for the shift operator T n : 



OO A 



(4) T n = ^ L - < = . 

i=0 

If u n depends simultaneously on the fast variable n and on the slow variable n\ i.e. u n = u(n, ni), 
the action of the total shift operator T n will give T n u n = ii n +i = u[n + l,n>i + e). So we 
can split it in terms of partial shift operators 7^ and 7^1 which are defined respectively by 
T n u(n,ni) = u{n + l,ni) and Tn\ u(n, n\) = u(n,n\ +e), T n = T n Tn) , where Tn) is also given 
by Eq. @. The dependence of u n on n can be easily extended to the case of one fast variable 
n and K slow variables rij = Ejn, Ej G 1, 1 < j < K, The total shift operator T n will now be 

written as: T n = 7^ ITfLi • 



Let us now consider a nonlinear partial difference equation 



(5) 



0, A/" (±) , 7W (±) > 0, 



for a function u n ^ m : 1? — > R which depends on two integer indexes n and m which we will call 
respectively space and £ime variables. Eq. ([5]) contains shifts of m contained in the intervals 
(m — m + A4^) and and n-shifts in the interval (n — Af^~\ n + AA + )). Under the 

hypotheses (a, b, c) we can give a series representation of the shifted values of u n ^ m around the 
point (n,m). Choosing 

e ni =JVie, e mj =Mje j , 1 < j < K, (s, Ni, Mj) e R 

we can write 

( e ) + °° N j 

(6a) T n = T n T ni 1 = Tn ^ ] ^ An ) >/ ^n = ~ T f^n 1 i 

3=0 °' 

K +oo 

(6b) T m = T m J] < mj = T m , 

j=l j=0 

K +oo 

f« r ^l tt — TT T^™ 1 ■* ITT j — TT r j A^ 

V U W ± n ± m — l n l m l n\ II 'mj — >n>m / c > Al n,m> 

3=1 j=0 

where the operators A m \ An}m are appropriate combinations of A m \ = -jrdm k and An) = 



jj-dni ( see |22| for more details). When we insert the explicit expressions ([6]) of the shift 
operators in term of the derivatives with respect to the slow variables into eq. ([5]) we get an 
e-dependent PDE of infinite order. Then we assume for the function u = u{n, m, ni, {rrij}j = i , e) 
a double expansion in harmonics and in the perturbative parameter e 

+oo 7 

(7) u{n,m,n 1 ,{m j }f =1 ,e)=J2 £ ( n l> Mf=i) K, m > 

7=1 a=— 7 
it. jl i{Kn-u{n)m] „(-«) — 

where the index 7 is chosen > 1 in order to let any nonlinear part of eq. ([5]) to enter as a 
perturbation in the multiscale expansion. 
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3.1.2. From derivatives to shifts. Splitting ([5]) in the various powers of e and in the different 
harmonics the multiple scale approach produces from a given partial difference equation partial 
differential equations for the amplitudes u^f 1 . Starting from the obtained partial differential 
equation we can write down a partial difference equation inverting the expression of the shift 
operator as 

(8) d ni = lnT ni =ln(l + h ni A«) = £ [ ~ >- Wl Affi', 

i=l 

where A^' = (T ni — l)/h ni is just the first forward difference operator with respect to the 
slow- variable n%. 

Only when we impose that the function u n is a slow-varying function of order I in n\, i.e. 
A^j 1 ^ 0, the d ni operator reduces to polynomials in A ni of order at most I. In [14], choosing 
/ = 2 for the indexes n\ and m\ and I = 1 for m,2, it was shown that the integrable lattice potential 
KdV equation [17] reduces to a completely discrete and local Nonlinear Schrodinger Equation 
(NLSE) which has been proved to be not integrable by singularity confinement and algebraic 
entropy |20,27J. Consequently, if one passes from derivatives to shifts, one ends up in general 
with a nonlocal partial difference equation in the slow variables n K and ms- 

3.2. The orders beyond the NLSE equation and the integrability conditions. The 

multiscale expansion of a difference equation ([5]) for analytic functions will give rise to a continuos 
PDE. So a multiple scale integrability test will require that an equation of the class of equations 
(PQ) is integrable if its multiscale expansion will go into the hierarchy of the NLSE. To be able to 
do so we need to consider here the orders beyond that at which one obtains for the harmonic 
the (integrable) NLSE. The first attempts to go beyond the NLSE order has been presented 
by Kodama and Mikhailov and by Santini, Degasperis and Manakov in (6j[T3]. Starting from S 
integrable models (models integrable via a Scattering Transform), through a combination of an 
asymptotic functional analysis and spectral methods, one succeeds in removing all the secular 
terms from the reduced equations order by order. The results could be summarized as follows: 

(1) The number of slow-time variables required to • • • for the amplitudes Ua s coincides 
with the number of nonvanishing coefficients C0j (k) = -p—^p-] 

(2) The amplitude evolves at the slow-times t a , a > 3 according to the a— th equation 
of the NLS hierarchy; 

(3) The amplitudes of the higher perturbations of the first harmonic j > 2, taking 
into account some asymptotic boundary conditions, evolve at the slow-times t a , a > 2 
according to certain linear, nonhomogeneous equations. 

Thus the cancellation at each stage of the perturbation process of all the secular terms is a 
sufficient condition to uniquely fix the evolution equations followed by every tt^ , j > 1 for 
each slow-time t a . Point 2 implies that a hierarchy of integrable equations always provide 
compatible evolutions for a unique function u depending on different times, i.e. the equations 
in its hierarchy are generalized symmetries of each other. In this way this procedure provides 
necessary and sufficient conditions to get secularity-free reduced equations [6]. 
Then, following Degasperis and Procesi [7] we state: 

Definition 3.1. A nonlinear PDE is said to be integrable if it possesses a nontrivial Lax pair 
and consequently an infinity of generalized symmetries. 

As a consequence of this Definition we have the following Proposition: 
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Proposition 3.1 If equation ([5|) is S integrable, then under a multiscale expansion the functions 
it,- , j > 1 satisfy the equations 



(9a) 
(9b) 

V j, a > 2, where K a 



a (!) 



,(!) 



M uf = U(j), M a = d to 



KL 



.(!) 



II 



(1) 



is the a-th flow in the nonlinear Schrddinger hierarchy. All the 



other Uj, k>2 are expressed in terms of differential monomials of Up \ p < j. 

In (|9b|) f a (j) is a nonhomogeneous nonlinear forcing term and K' a [u] v is the Frechet derivative of 
the nonlinear term if CT [it] along the direction v K' a [u]v = j^K a [u-\-sv] \ s =o, i- e. the linearization 
of K a [u] along the direction v. Eqs. ([9]) are a necessary condition for integrability and represent 

a hierarchy of compatible evolutions for the same function at different slow-times. The 
compatibility of eqs. (l9b|) is not always guaranteed but is subject to a sort of commutativity 
conditions among their r. h. s. terms f a (j). Then it is easy to prove that the operators Mo- 
defined in eq. (|9bj) commute among themselves. Once we fix the index j > 2 in the set of eqs. 
(|9bp . this commutativity condition implies the following compatibility conditions 

(10) M <T f a ,(j)=M a ,f a (j), Va,a'>2, 

where, as f a (j) and f a i (j) are functions of the different perturbations of the fundamental har- 
monic up to degree j — 1, the time derivatives dt a , df , of those harmonics appearing respectively 
in M a and M a i have to be eliminated using the evolution equations Q up to the index j — 1. 
The commutativity conditions f 1 1 [) turn out to be an integrability test. 

Following [6] we conjecture that the relations ([9]) are a sufficient condition for the integrability 
or that the integrability is a necessary condition to have a multiscale expansion where eqs. ([9j) 
are satisfied. To construct the functions f a (j) according to the Proposition 3.1 we define: 



Definition 3.2. A differential monomial p 



(i) 



j > 1 in the functions it,- , its complex 



conjugate and its ^-derivatives is a monomial of "gauge" 1 if it possesses the transformation 
property 



P 



e i9 p 



ii 



(i) 



r(i) 



e ie uf\ 



Definition 3.3. A finite dimensional vector space V n , n > 2 is the set of all differential poly- 
nomials in the functions uf \ j > 1, their complex conjugates and their ^-derivatives of order n 
in e and gauge 1 when 



order [d'fUj 



(i) 



order ( d. 



'£ ll Uj ) = m + 3i > ri - (): 



Definition 3.4. V n {yn), m > 1 and n > 2 is the subspace ofV n whose elements are differential 
polynomials in the functions uj 1 ' 1 s, their complex conjugates and their ^-derivatives of order n 
in e and gauge 1 for 1 < j < m. 

From the definition (|3.4p one can see that in general K a 



u 



(i) 



G d^u^ U 7Vk(1) and that 
fcr(j) € Vj+a-(j — 1) where j, a > 2. The basis monomials of the spaces V n (m) can be found, 
for example, in [22] 

Proposition 3.2 If for each fixed j > 2 the equation M0\) with a = 2 and a' = 3, namely 
■^2/3 (j) = M3/2 (j), is satisfied, then there exist unique differential polynomials f a (j) Vcr > 4 
such that the flows M a Uj = f a (j) commute for any a > 2. 
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Hence among the relations (jlOp only those with a = 2 and a' = 2 have to be tested. 
Proposition 3.3 The homogeneous equation M a u = has no solution u in the vector space 
V m , i.e. Ker{M a )C\V m = 0. 

Consequently the multiscale expansion ([9]) is secularity-free. Finally we define the degree of 
integrability of a given equation: 

Definition 3.5. If the relations U0\) are satisfied up to the index j, j > 2, we say that our 
equation is asymptotically integrable of degree j or Aj integrable. 

3.2.1. Integrability conditions for the NLSE hierarchy. We specify here the conditions for as- 
ymptotic integrability of order k or integrability conditions. To simplify the notation, we 
will use for Uj the concise form u(j). Moreover, for convenience of the reader, we list the fluxes 
K a [u] of the NLSE hierarchy for u up to a = 4: 



(11a) 
(lib) 

(11c) 

(lid) K 4 [u] = -\C I Utftf + 



K 2 [u] 



K\ [u] = Au£ 
-ipi 



, P 2 \ |2 

Ucc H \u\ u 

pi 



KM 

P2 



B 



U£££ H \U\ Uc 

Pi 



Pi 



3y^?2 i|4 ii2 2— ii2 2 — 

\u\ u + 4\u\ nee + Sucu + 2\ut\ u + u net 

2pi 



where pi, P2, A, B and C are arbitrary complex constants. 

The A\ integrability condition is given by the reality of the coefficient p 2 of the nonlinear 
term in the NLSE. 

The A2 integrability conditions are obtained choosing j = 2 in the compatibility conditions 
([TP]) with a = 2 and a' = 3 

M 2 h U) = M 3 f 2 (j) . 



(12) 



In this case we have that / 2 (2) G P 4 (l) and / 3 (2) G V 5 (l) with dim(7 :, 4(l)) = 2 and dim(7 :, 5(l)) = 
5, so that f 2 {2) and /a(2) will be respectively identified by 2 and 5 complex constants 

(13a) f 2 (2) = cm e (l)Kl)! 2 + 6n ? (l)n(l) 2 , 

(13b) h{2) = a\u{l)\ 4 u{l) + /3|u e (l)| 2 u(l) + 7^(l) 2 n(l) + 5%(l)n(l) 2 + e|n(l)| 2 u^(l). 

In this way, if /?2 7^ 0, eliminating from eq. (|12p the derivatives of u(l) with respect to the 
slow-times t 2 and *3 using the evolutions (f9a|) with <r = 2 and cr' = 3 and equating term by term, 
we obtain the A 2 integrability conditions 



(14) 



a = a, b = b. 



So at this stage we have two conditions obtained requiring the reality of the coefficients a and 
b. The expression of a, /3, a, 5 in terms of a and 6 are: 



(15) 



a 



3iBp 2 a 







3LB6 
Pi 



7 



3i-Ba 



0, 



7- 



The ^3 integrability conditions are derived in a similar way setting j = 3 in eq. (I12p . In this 
case we have that / 2 (3) G P 5 (2) and / 3 (3) G V 6 (2) with dim(P 5 (2)) = 12 and dim(P 6 (2)) = 26, 
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so that /2(3) and / 3 (3) will be respectively identified by 12 and 26 complex constants 

/ 2 (3) = n|n(l)| 4 U (l) + r 2 |u ? (l)| 2 n(l) + r 3 |n(l)| 2 n^(l) + r 4 %(l)u(l) 2 + r 5 n f (l) 2 u(l) + 
+t 6 ^(2)|u(1)| 2 + r 7 u 5 (2)n(l) 2 + r 8 u(2) 2 u(l) + r 9 |n(2)| 2 n(l) + ri «(2)^(l)ii(l) + 

(16a) +rii«(2)fi € (l)u(l) + r 12 n(2)^(l)u(l), 

/ 3 (3) = 7l |n(l)| 4 ^(l) + 72 | U (l)| 2 n(l) 2 n f (l) + 73 | U (l)| 2 n^(l) + 74 u(l) 2 % f (l) + 

+75|^(l)| 2 «e(l) +76%(l)^(l)^(l) +77^(l)%(l)n(l) +7 8 u ?? (l)u 5 (l)u(l) + 
+ 79 |u(l)| 4 n(2) + 7 io|u(l)| 2 n(l) 2 u(2) + 7 n%(l)«(2) 2 + 7 i2^(l)|u(2)| 2 + 
+ 7l3 |^(l)| 2 u(2) + 7l4 |n(2)| 2 u(2) + 7l5 n ? (l) 2 u(2) + 7l6 |u(l)| 2 ^(2) + 
+ 7 i 7 u(l) 2 ^(2) +7i a ti(2)« tt (l)«(l) + 7 i9«(2)ti K (l)«(l) + 720 u(2)^(l)u(l) + 
^iti(2)u € (2)«(l)+Tto«(2)u C (2)u(l)^ + 

(16b) + 725 %(2)« e (l)«(l) + 726 ^(2)u(2)«(l). 

Let us eliminate from eq. f)12|) with j = 3 the derivatives of u(l) with respect to the slow-times 
t 2 and i 3 using the evolutions (j9a|) respectively with a = 2 and <r' = 3 and the same derivatives 
of u{2) using the evolutions ()9b|) with cr = 2 and a' = 3. Let us equate the remaining terms term 
by term, if p 2 ^ 0, and indicating with E4 and Ii the real and imaginary parts of Tj, i = 1, . . . , 12, 
we obtain the A 3 integrability conditions 

aJ 6 (6 - a)J 6 al 12 R 2 (a- b)I 6 al 12 

4pi 2,o 2 2p 2 2 4p 2 4p 2 

fl 2 , (o-b)/ 6 , (2b-o)Ji2 a/ 8 (a-b)J 8 

-K5 — — H 1 , Kq — , tii — -K12 H 5 

2 4p 2 4p 2 p 2 p 2 

Rs = R9 = 0, -Rio = -R12, = -R12 H ~, 

P2 

_ (6 + a)fli 2 , pi/i , I 2 -h-2h , [26(a - 6) + a 2 ] J 8 

J 4 = : 1 1 : 1 7-9 ; h — U, 

4p 2 p 2 4 4^ 

(17) / 9 = 2/ 8 , I W = I 12 , hl = h + h2- 

For completeness we give the expressions of the jj, j = 1, . . . , 26 as functions of the Tj, i = 
1,...,12: 



71 



72 



3B 
8pJ 
3Ba 



2a 2 I 

2bR 12 - 8 Pl h + 2(/ 2 - 2/3 - 2I 5 )p 2 + i(6 - 5a)/ 6 + - 3ia/ 12 

P2 



2 



4p( 



(a - 2b)h _ ■ 

1^6 H 1" 7"12 

P2 



3ii?T4 3i-Br5 3B(p 2 Iq + 3ai/ 8 ) 

76 = , 77 = 75 , 78 = 73 + 





3iBr 3 


, 73 


2pi ' 


3iBr 5 






79 = - 


Pi 




3iBr 9 


3 


2pi ' 


713 = 



3iBr 2 

7 4 = 0, 75 



2p 



1 

- 3ai/ 8 ) 

pi pi M 

3iBp 2 R d 3iBr 9 3iBr n 3\Bt 12 

710 = — 7T2 — ' 7n=0, 7 i2 = — , 713 = —^ , 714 = 0, 715 = — . 

2p{ 2pi 2pi 2p 1 

3iBtq 3iBtiq 3iBt$ 
716 = -7, , 717 = 718 = 0, 719 = —?; , 720 = 7i5, 721 



2pi ' ' 2pi " " pi 

3iBri 

(18) 7 22 = 712, 723 = 716 + 719, 724 = 713, 725 = , 726 = 0. 

pi 

The conditions in the case of C integrable or linearizable equations are similar to those 
presented here and can be found in ref. |23j . 
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4. Multiple scale Expansions of the exceptional models 

Let us perform the multiscale analysis of the differential-difference models ([1]) for the real 
function u n (t) with the functions Q n given by ([2]). To do so we expand the variable u n according 
to the prescriptions of the previous section and split the obtained system in terms of the different 
harmonics and various powers of e. The equations at the various orders are presented in the 
Appendix A, here we just present the final results. 

Proposition 4.1 All models ^ pass the A± and A2 integrability conditions however fail the A3 
one for any choice of the involved parameters ^ . . 

Moreover, even the particular choice done by Barashenkov, Oxtoby and Pelinovsky in the 
case of Qi and Q2 is no more integrable than the others. 

In Appendix A it is also shown that these nonlinear model can never be nontrivially linearized. 

5. Concluding remarks 

In this article we have shown, applying the multiple scale expansion integrability test, that the 
set of discrete 4 models constructed by Barashenkov, Oxtoby and Pelinovsky are not integrable 
but have just some nontrivial degree of integrability as they pass the two lowest integrability 
conditions. 

The models do not satisfy the symmetries of the continuous </> 4 theory but just possess 
travelling kink solutions. It would be interesting to construct discrete models which preserve the 
symmetries of the continuous one. This can be done [16], however this implies that we will need 
to have a variable non uniform lattice. As the continuous model is not completely integrable, 
we do not expect such discrete symmetry preserving model to be completely integrable. 

An interesting extension of the calculation done here consists in expressing the multiple scale 
expansion in term of normal forms. In such a way we would be able to use the results to give 
approximations to the solutions of the starting equations. Work on both these aspect is in 
progress. 
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Appendix A. Details of the calculations 

We present here in the following the equations one obtains by splitting (pQ) , once we substitute 
the field u n by its multiple scale expansion ([7]), into the various harmonics and different orders 
of e. Except when explicitly stated, the equations are valid for any function Q n given in ([2]). 

• Order e and a = 0: One obtains 

(19) U <°> = 0; 

• Order e and a = 1: If one requires that 7^ 0, one obtains the dispersion relation 



, _ 4sin 2 (^/2) 1 
(20) u p -. 

• Order e 1 and a = 0: We obtain 



(21) 



4 0) = 0; 
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(22) d tl u?-a ldni u? = 0, ai ^ sin ^) 



h 2 



which tell us that has the form 

, \ fil \ * . , sin (kH) 

(23) u\ L) =g(Z,t j ,j >2), g = /mi + A- ^ 

where g is an arbitrary function of its arguments going to zero as £ — > ±oo; 

• Order e 2 and a = 2: Taking into account the dispersion relation (|20p . we have 

(24) 4 2) = 0. 

• Order e 3 and a = 0: Taking into account relation (fl~9|) . we have 

(25) 4 0) = 0; 

• Order e 3 and a = 1: Taking into account the dispersion relation ([20]) and the eqs. (fl~9 
[22]) . we have 



(26a) d tl - aid ni vtp = -d t2 u^ - ipid^u^ - ip 2 \u^\ 2 , 

3 + (h 2 - 4) cos (nh) + cos (2/e/i) 
( ^ Pl = 2^ [/i 2 -4 + 4cos(/c/i)] ' 

. . 4^2 + (5ui + 3/X3) cos (k/i) + 2/i2 cos (2k7i) + u\ cos (3k/i) 
(26C) P2 = 20^ ' 



(26d> 2 



(26e>J2 



(5/i 2 + 8/U3) Hi - + [(11^1 - 5 ^2) Ml + 3 (A*2 + ^3) ^3] cos (nh) 



CO 



+ 4//i/i3 - 2/i 2 ^3) cos (2k/i) + /Ui (izi - ti 2 ) cos (3k/i) 



Co' 



3 (4 + h 2 ) p 1 + 2/i 2 + 3 (1 - 2/x 2 ) cos (nti) + [3 (4 + h 2 ) m + 4/i 2 ] cos (2/cfe) 



4w 



where (|26cP is obtained in the case of the model Q lf ([26d]) for Q 2 and l|SSe|l for Q 3 . As 
a consequence of eq. (|22p . the right hand side of eq. (|26ap is secular. Hence we have to 
require that 



(27a) 
(27b) 



^ 2 



(i) 



Eq. (|27a|) tells us that u% also depends on £ while (j27bj) is an integrable NLSE, as 
from the definitions (I26cl I26d[ I26ef) we can see that p 2 is a real coefficient in all of the 
three cases. Hence all the models are Ai— integrable. 
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If we require that our models be A±— linearisable [23], then we need p 2 = for any k. 
In this case we have respectively 



(28a) 
(28b) 
(28c) 



Qi 

Q2 



Ml = A*2 = M3 
fj, t = fj, 3 = 0, 
no solution, 



0, 



i.e. in the first two models only the trivial linear cases are selected while the third model 
doesn't admit an A\— linearisable reduction. 

Order e 3 and a = 2: Taking into account the dispersion relation ([20]) and (fT9l |24"|) . we 
have 



(29) 



4 2) = o; 



Order e 3 and a = 3: Taking into account the dispersion relation ([20]) . we obtain 



(30a) 
(30b) 

(30c) 
(30d) 



.(3) 



(1)3 

a 2 u\ , 



«2 



a 2 



a 2 



. /i 2 [fj,^ + 2/Lt2 cos (k/i) + 2/^i cos (2k/i)] cos (nh) 
~ 20 [8 - 2h 2 - 9 cos (re/i) + cos (3kR)] ' 
^_ h? [fi 3 + 2/ii cos (ah)] [fii - fi 2 + (a*2 + M3) cos + y,\ cos (2k/i)] 
8 — 2h 2 — 9 cos (nh) + cos (3nh) 
h 2 [(4 + /i 2 ) /ii + 2/i 2 + (1 - 2/U 2 ) cos (k/i) + (4 + h 2 ) m cos (2k/i)] 
4 [8 - 2/i 2 - 9 cos (k/i) + cos (SkH)] 



where (|30bp is obtained in the case of the model Q\, ()30cp for Q 2 and (|30d|) for Q3. 
Order e 4 and a = 0: Taking into account (fl~9l I2TI [24"|) . we get 



(31) 



,(°) 



0; 



Order e 4 and a = 1: Taking into account the dispersion relation (|20p and eqs. (fT9 
IZ7a |l57E|> . we get 



(32a) 
(32b) 



« „(i) fl „(i) 



5t 3 «i 



(i) 



.(1) 



(i) 



/ 2 (2; 



D 



[6 - 8/1 2 + h 4 + 2 (/i 2 - 4) cos (k/i) + 2 cos {2nh)] sin (k/i) 



6/i 2 w [/i 2 -4 + 4 cos (kTi)] 



(1) 



(1) 



-1P1 



9 | n W + £» ( U « 2 4D + 2 \u\ 



Pi 



In the above relations 



(i) 



(jllcp is the cmKdV flux, the second flux of the NLSE 



hierarchy, K' 2 



,(!) 



is the Frechet derivative along the direction of the NLSE 
defined by relation ()27b|) and / 2 (2) (|13a|) is a nonlinear forcing term 
depending on and defined by the coefficients a and 6. As a consequence of (|22l I27aj) 



flux if 2 



(i) 



(33a) 
(33b) 



(34a) 
(34b) 

(35a) 



(35b) 



(35c) 



(35d) 
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the right hand side of (|32ap is secular, so that 



ii 



" ^2 



ai9 ni «^ = 0, , 



,(!) 



+ /a(2). 



Eq. ()33ap tells us that depends on £ too while in (|33b|) . as a consequence of (|27b|) . 
the first term of the right hand side is secular, so that 



d t3 u^ - K 3 



4 1) = / 2 (2) 

= 0. 



The coefficients a and b of the forcing term fi (2) (|13a|) are given by 

[-36 (ft 2 - 4) mi + 3 (4 - 8ft 2 + ft 4 ) M2 - 5 (ft 2 - 4) ^t 3 ] sin (reft.) 

a= a; + 

[(-153 + 32/i 2 - 4ft 4 ) Ml - 6 (ft 2 - 4) /j 2 - (2 + 8ft 2 - /i 4 ) Ms] sin (2reft) 

+ 

[23 (h 2 - 4) jiii + (36 - 8ft 2 + /i 4 ) M2 - 3 (ft 2 - 4) 113] sin (3reft) 
A a 

[(-42 + 8ft 2 - ft 4 ) m - 3 (ft 2 - 4) M2 + 5^ 3 ] sin (4reft) 

+ aI 

3 (ft 2 — 4) /ii sin (5reft) 4- Ml sin (6reft) 
A^ ' 

A a —40a; ft —4 + 4 cos («ii)J [3 + (ft 2 - 4) cos (reft) + cos (2reft)l , 

_ {l2/i 2 + (ft 2 - 4) M3 + 2 [7/ii + (ft 2 - 4) 112 + 5^3] cos (reft) 
~ At 

2 [(ft 2 — 4) Ml + 4/i2] cos (2reft) + 6mi cos (3reft)} sin (reft) 
A5 ' 

A b =20u [ft 2 - 4 + 4 cos (reft)] , 

{ [(50 - 48ft 2 + 7ft 4 ) ii, 2 - 46 (ft 2 - 4) Mi + 6 (4 - 8ft 2 + ft 4 ) M.-i Mi + 

a = 

A a 

+ [(70 - 21ft 2 + 2ft 4 ) H2 - 5 (ft 2 - 4) /x 3 ] M3} sin (reft) 

A^ + 
{ [(1 + 10ft 2 ) M2 + (-157 + 16ft 2 - 2ft 4 ) mi] m + 
+ A^ 

+ [(-58 + 6ft 2 + ft 4 ) 112 - 12 (ft 2 - 4) m + (-2 - 8ft 2 + ft 4 ) 113] M3} sin (2reft) 

A a 

{- [(59 - 19ft 2 + ft 4 ) 112 + 17 (ft 2 - 4) tlx +2 (36 - 8ft 2 + ft 4 ) 113] Mi + 
+ Al 

+3 [(ft 2 + 4) ii 2 + (ft 2 - 4) 113] fis} sin(3reft) 

Al + 
{[(46 - lift 2 + ft 4 ) fi 2 - (32 - 8ft 2 + ft 4 ) in - 6 (ft 2 - 4) 113] m + 5 (M2 + M3) M3} sin (4reft) 

+ aI 

{ [(13 — 3ft 2 ) 112 + 3 (ft 2 — 4) Ml] Ml + 2M2M3} sin (5reft) + (mi — M2) Ml sin (6reft) 

Al ' 

A a =2w [ft 2 - 4 + 4 cos (reft.)] [3 + (ft 2 - 4) cos (reft) + cos (2reft)j , 

{ [(ft 2 - 4) in - 3 (ft 2 - 7) IJ.2 + 16^3] Ml + [(ft 2 - 8) IJ.2 + (ft 2 - 4) ii 3 ] 113} sin (refc) | 



b = - 



{2 (7mi - 6M2)Ml + [2 (ft 2 - 4) Ml + (13 - 2ft 2 ) M2 + 5 113] 113} sin (2reft) 
+ At. + 

{ [(ft 2 - 4) Ml - (ft 2 — 5) M2] Ml + 2 (4mi - 3M2) 113} sin (3reft) + 3 (mi — M2) Ml sin (4reft) 



A„ 



A b =ui [ft 2 - 4 + 4 cos (reft.)] 
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{4 (-5 + 58/ii + 35p 2 ) + h 2 [5+ (26 - 4ft 2 + ft 4 ) pi + 2 (2ft 2 - 21) fi 2 ] } sin (reft) 

a — — 

A„ 

[-2 + 352jUi + ft 4 (1 - 22/ji - 2fi 2 ) + 116^2 - 4ft 2 (2 + 3/i 2 )] sin (2reft) 
+ Al 

{12 (1 + 19/ii + 2p 2 ) + ft 2 [-3 + (25 - 4ft 2 + /i 4 ) m + 6/x 2 ] } sin (3reft) 

Al + 

[5 - 3 (ft 4 - 16) (Hi - 10/j 2 ] sin (4reft) + [(ft 2 + 4) fii + 4^ 2 ] sin (5reft) 
+ Al ' 

A a =8w [ft 2 - 4 + 4 cos (reft)] [3 + (ft 2 - 4) cos (nh) + cos (2reft)] , 

_ [ft 2 (1 + 7/n - 2^2) + 4 (-1 + 7/ii + 4^2)] sin (reft) 
~ A5 

[5 + 2 (ft 4 - 16) Mi + 2 (2ft 2 - 13) fj, 2 ] sin (2reft) + [7 (ft 2 + 4) fi X + 12/i 2 ] sin (3reft) 



A 6 =4w [ft 2 - 4 + 4 cos (reft)] , 

where (|35al I35b[) is obtained in the case of the model Qi, (|35cl I35d|) for Q2 and (|35el 
35f[) for Q3. As one can see, in all the three cases the coefficients a and b are real. As a 
consequence all the three models are A2 — integrable; 

Order e 4 and a = 2: Taking into account the dispersion relation (|20p and (|19[ [2T| [241 



29J), it results 



4 2) = 0; 



Order e A and a = 3: Taking into account the dispersion relation (I20p and (|19[ l22l [30]) . 
it results 



u 



(3) 



a 3 U2 + 



( 4 A J 1 ) 2 . 



«J I « 



'1 ! 



(3) 

As do not enter into the final result, the coefficients 03 and 04 are not explicitly 
written down here. 

Order e 4 and a = 4: Taking into account the dispersion relation (|20p and (|24p . it results 

4 4) = 0. 

Order e 4 and a = 0: Taking into account ([19l I2T1 [2il 1251 1291) . we get 



(o) 



Order e 5 and a = 1: Taking into account the dispersion relation (|20p and (|19| [2T| [22l 
Ml EZH [270 SJ3 [33al IMallMbj) . we get 



d tl u¥ - a ldni u? 



X 1 ) 



,(!) 



/ 2 (3) 



t 3 "2 



/s(2; 



u 



(1) 



^ 7 (5 - 8h 2 + /i 4 ) + (h 2 - 4) (14 - 8h 2 + /i 4 ) cos (re/i) 
2Ah 2 uj[h 2 -4 + 4cos(k/i)] 3 
(28 + 8/1 2 - /i 4 ) cos (2k/i) + 2 (/i 2 - 4) cos (3k7») + cos (4/sfc) 
24/i 2 w[/i 2 -4 + 4cos(k/i)] 3 



In the above relations 



?/ 



(i) 



(jlldp is the is the third flux of the NLSE hierar- 



chy, K> 



.(!) 



U2 is the Prechet derivative along the direction of the cmKdV flux 
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v 



(i) 



defined by relation (jllcp and fi (3), /3 (2) are the nonlinear forcing terms de- 
fined in (|13b| I13aj) . As a consequence of ([22l I27al I33ap the right hand side of (|40ap is 
secular, so that 



(41a) 
(41b) 



d t2 u 3 



(i) 



K> 2 



.(!) 



'*4 U 1 



(1) 



K' 3 
K 4 



(1) 



f ' " /s (2) 
+ / 2 (3). 



The first relation tells us that u^p depends on £ too while in (|41bp . as a consequence of 
dMa|) and of 



.(1) 



/ 3 (2)= 5 t3 -^ 



/2 (2) , 



the first term on the right hand side is secular. Moreover also the second term on the 
right hand side is secular because, when we equal it to zero, we obtain a generalized 
symmetry the NLSE, the third equation of the corresponding hierarchy. Hence 



(42a) 
(42b) 



d t2 uP - K' 2 
d t3 uP-K' 3 



u 



(1) 



,(!) 



,(!) 



/ 2 (3), 



d u uP - K 4 



,(!) 



+ /s(2). 



Finally the first term on the right hand side of (|42bp is secular because, when we equal 
it to zero, we obtain a generalized symmetry the cmKdV equation as both equations 
belong to the same NLSE hierarchy. Hence 



(43a) 
(43b) 



d t4 u { p - K A 



u 



(i) 



u 



(1) 



4 1} = h (2) , 
= 0. 



The real and imaginary parts (Rj, Ij) of the coefficients Tj, j = 1,. . . ,12 of the forcing 
term /2(3) are given by 



(44) 
(45) 



Ri — R>2 — R3 — Ra — R5 — Rs — R9 — 0, -Rio — R12, 
h = I7 = ho = hi = I12 = 0, Iq = 2/g, 
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independently from the model. 



(46a) Ji 

(46b) I 2 = 

(46c) I 3 = 

(46d) / 4 = 

(46c) J 5 = 

(46f) Re = 

(46g) fi 7 
(46h) J 8 

(46i) fl 10 

(46j) flu 



(47a) 



(47b) 



2cs 2 M2Pi — 5 (2pjp2 + 3aBw — 6Cp2") p 2 + a 2 Pi {p2 [3 cos (2nh) + cos (4rc/i)] + 

20p\oj 

+ [p3 + 2pi cos (2kK)\ [2 cos (k/i) + cos (3t^h)\ } 



20pfw 




2p 2 pi + 20 [(a + 26) aipi + (3Bai 


- 2p\ + 2Cu>) , 


20piw 




+6pi cos (k/i) — p2 cos (2nh) — 2pi 


cos (3k/i) 


20ui 




2aaipi + 2 (3Bai - 2pJ + ACoj) pi 


— 3aBco-\- 



2p\oj 

+ [4pi + P3 + 3p2 cos (k/i) + 2pi cos (2k/i)] cos (/t/i) 
20w 

2 [(206ai + JU2) pi + 20Cp2^] + [(5jUi + P3) cos (k/i) + P2 cos (2nh) + pi cos (3k/i)] pi 

40piw 

20 (2aipi - 3Boj) a - p 2 pi + 40 (3Bai - p\ + 3Cu>) p 2 + 2 [p 2 + 2m cos (reh)] pi cos (2k/i) 

40piw 

(2ctipi — 3Bw) P2 [4pi + P3 + 3p2 cos (kH) + 2pi cos (2nh)\ sin (k/i) 

Pilj 10w 
20aip2 + [p3 + 2pi cos {2kK)\ sin (k/i) 4- P2 sin {2kK) 
20bj 

(5pi + 3^3) cos (nil) + 2 [2 + cos (2kH)] P2 + Pi cos (3«/i) 
20w 

(2aipi — 3Buj) p2 [4pi + P3 + 3p2 cos (k}i) + 2pi cos (2nh)\ sin (k/i) 

Pilj IOcj 
20aip2 + [P3 + 2pi cos (2nh)] sin (k/i) + P2 sin (2/c/l) 



for the model Qi 

4 (3pip 2 + 4pip 3 - P2M3) c<2P? - (2p?p2 + 3aBu - 6Cp 2 aj) p 2 + 

7l = 7~2 

4pjw 

+4a 2 pi {[(7pi - 2p 2 ) Pi + 2 (p 2 + p 3 ) P3] cos (nh) + 3pi (p 2 + 2p 3 ) cos (2nh) + 

ipjui 

+ [(4pi - 3p 2 ) pi + (P2 + P3) P3] cos (3nh) + (pi - p 2 ) [2p 3 cos (4k h) + pi cos (5/tfe)]} 

4pfw 

(a + 26) aipi + (3Bai - 2p\ + 2Cu) p 2 - 36Bw + 2 (p 2 + 2p 3 ) P1P1 + 

12 = 

piw 

+2 {3pi cos (nil) + (p2 — Pi) [P3 cos (2Kh) + pi cos (3kH)}} 



(47c) I 3 = 



(47d) I 4 



2aaipi + 2 (3Sai - 2/>? + 4Cw) p 2 - 3aBu+ 
2piai 

+ {2 (3pi - P2) Pi + (P2 + P3) P3 + 2 [(p 2 + 3p 3 ) pi - P2P3] cos (kK) + 2 (pi - p 2 ) Pi cos (2nh)} cos (nh) 

ui 

2 [6«i + (P2 + 2p 3 ) pi] pi + 2Cp 2 ^ + {[(7pi - 3p 2 ) pi + (P2 + P3) P3] cos (nh) + 

2piaj 

+ (pi — P2) [2p3 cos (2kK) -f- Pi cos (3^6)]} pi 



(47e) h 



(47f) R 6 = 



2p\LO 

6 {Bai + Clu) p 2 - 2 [(pi - p 2 ) P3 + P1P2] pi + (2ceipi - 3Bu) a+ 
2piw 

+2 [(pi + 2p 2 ) cos (nh) + (p 2 + 2p 3 ) cos (2nh) + (pi - p 2 ) cos (3ft/i)] pipi 

2piai 

(2aipi - 3Buj) p 2 - 2 {2 (3pi - p 2 ) Pi + (p2 + P3) P3 + 
piw 

+2 [(p 2 + 3p3) pi — P2P3] cos (nh) + 2pi (pi - p 2 ) cos (2nh)} pi sin (nh) 
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Q1P2 + [(pi - 3p 2 ) Ml + (M2 + M3) M3] sin (k/i) + (pa - p 2 ) [2p 3 sin (2«fc) + Hi sin (3k.Ii)] 
(48a) R 7 = 1 



(48b) 



(49a) 



(49b) 



. (5p 2 + 8/13) AH — M2M3 + — 5^2) Ml + 3 (M2 + M3) M3] cos (k/i) + 

— 

u 

+ [(jl2 + 4p 3 ) hi - 2p 2 p 3 ] cos (2k/i) + (/ii — p 2 ) Ml cos (3k/i) 



, .„ s „ (2aipi - 3-Bw) P2 - 2 {2 (3pi — M2) Ml + (M2 + M3) M3 + 

(48c) i?i = 

piu> 

+2 [(M2 + 3M3) Ml — M2M3] cos (k/i) + 2 (pi — M2) Ml cos (2k/i)} pi sin (k/i) 

piU 

2 {aip 2 + [(Ml — 3M2) Ml + (M2 + M3) M3] sin (k/i) + (fii - H2) [2p3 sin (2k/i) + Ml sin (3k/i)]} 
(48d) Rn = . 

for the model Q2 and 

[(/i 2 + 4) Hi + 2 M2] o 2 pl - (2p 2 p2 + 3aBui - 6Cp2u) P2 + «2P? {2 (1 - 2p 2 ) cos (k/i) + 

^1 — ~7~2 

4pjw 

+3 (h 2 + 4) pi cos (2k/i) + (1 - 2p 2 ) cos (3k/i) + [(/i 2 + 4) pi + 2p 2 ] cos (AkK)} 

ipjui 

_ 2 [(a + 2b) aipi + (3Bqi - 2p\ + 2Cw) p 2 - 36Sw] + 



(49c) I 3 = 

(49d) h = 

(49c) / 5 = 

(49f) R 6 = 

(49g) fl r = 

(49h) I 8 = 

(49i) Rid 







2piLj 


+ (ft 2 H 


- 4) hi - 


4piP2 - [(/i 2 + 4) pi + 2p 2 ] cos (2k/i) 






2lu 


2ao:ipi 


+ 2 (35 


ai - 2p 2 + 4Cu) p 2 - 3aBw+ 






2piu 




2p2 + 4 


[{h 2 + 4) pi + p 2 ] cos (k/i)} cos (k/i) 






4w 


2 [4bai 


+ (/! 2 + 


4) pi] pi + 8Cp 2 w + { (1 - 2p 2 ) cos (k/i) + 2 



8piu 

12 (Bai + Cuj) p 2 - [(/i 2 + 4) pi + 2 (p 2 + 2pip 2 )] pi + 2 (2aipi - 3Buj) a + (h 2 + 4) pipi cos (2k/i) 

4piw 

(2aipi - 3Bui) p 2 (2p 2 - 1) sin (k/i) - 2 [(h 2 + 4) pi + p 2 ] sin (2k/i) 



(49j) R„ 



Pilo 2lu 
4ctip 2 + (1 - 2p 2 )sin (k/i) + 2 [(h 2 + 4) pi + 2p 2 ] sin (2k/() 
4ui 

3 (h 2 + 4) pi + 2p 2 +3(1- 2p 2 ) cos (k/i) + [3 (h 2 + 4) pi + 4p 2 ] cos (2k/i) 

4a> 

(2aipi - 3Bic) p 2 (1 - 2p 2 ) sin (k/i) + 2 [(/i 2 + 4) pi + M2] sin (2k/i) 

Pilo 2u 
4aip 2 + sin(K/i)[l - 2p 2 ] + 2sin(2K/i)[(/i 2 + 4)pi + 2p 2 ] 



for the model Q3. As the coefficients (pl6l WI\ H5 1 H9|) do not satisfy the algebraic relations 
(fT7|) . for any Q the A3 integrability is never satisfied. 
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